A Kac-potential treatment of nonintegrable interactions 
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We consider d-dimensional systems with nonintegrable, algebraically decaying pairwise interac- 
tions. It is shown that, upon introduction of periodic boundary conditions and a long-distance cutoff 
in the interaction range, the bulk thermodynamics can be obtained rigorously by means of a Kac- 
potential treatment, leading to an exact, mean-field-like theory. This explains various numerical 
results recently obtained for finite systems in the context of "nonextensive thermodynamics," and 
in passing exposes a strong regulator dependence not discussed in these studies. Our findings imply 
that, contrary to some claims, Boltzmann-Gibbs statistics are sufficient for a standard description 
of this class of nonintegrable interactions. 
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I. INTRODUCTION 



In studies of critical phenomena the range of the pair- 
wise interaction that couples the degrees of freedom is 
an important consideration. For interactions which de- 
cay algebraically at large distances, three classes of crit- 
ical behavior may be obtained. With the standard no- 
tation u(r) ~ — l/r d+a , one finds for system dimension- 
ality d < 4 that the criticality may be characterized as 
short-range for a > 2 — rj SI , non-classically long-range 
for d/2 < a < 2 — r] SI , and classically long-range for 
< a < d/2, where ?7 sr is the correlation-function ex- 
ponent in the corresponding system with short-range in- 
teractions (3,0] ■ The critical behavior matches at the 
bordering cases (e.g., a — ► d/2 from above and below) 
with additional logarithms. However, for a < the in- 
teractions are nonintegrable, i.e. Jd d ru(r) — ► oo, and 
so, under standard definitions, the thermodynamic limit 
does not exist. (See Refs. J3|-|| for rigorous treatments.) 

Nevertheless, recent studies have focused on this non- 
integrable regime, typically using a finite system size to 
render the total system energy finite. These results are 
then interpreted as "nonextensive thermodynamics" QS|, 
in which the system energy density scales with some 
positive power of the system size, as do intensive vari- 
ables such as temperature. Examples of such works in- 
clude molecular-dynamics simulations of two- and three- 
dimensional systems with variants of the Lennard- Jones 
potential |7]|| , Monte Carlo simulations of one- and two- 
dimensional Ising PJlQ] and Potts Jll| systems, and a 
numerical study of the XY chain jl^ . Monte Carlo sim- 
ulations have also suggested classical critical behavior 
for a stochastic cellular automaton with long-range in- 



teractions in the regime a < Jl3| |. On the basis of 
these numerical studies of finite systems, several authors 
have conjectured that "nonextensive criticality" should 
be classical [p 11 12[|. 

Here we present an alternate approach to noninte- 
grable interactions, quite similar to, and in certain cases 
equivalent to, the use of the Kac potential. We introduce 
a long-distance cutoff in the interactions at some finite 
distance R, which enables us to use periodic boundary 
conditions and thus consider homogeneous systems. We 
then find the energy density in the thermodynamic limit 
to scale as a power of R rather than of the system size L, 
that is, we maintain extensitivity. By multiplying the 
pair interaction by the appropriate negative power of R 
we recover a well-defined R — > oo limit. In this way, 
by using standard methods (including the conventional 
canonical ensemble), we find an exact solution for the 
free energy for all — d < a < 0, and so demonstrate the 
classical nature of the criticality without resort to simu- 
lations or conjecture! These results have been announced 
previously in Ref. E3| . 

Indeed, this is not surprising — the "infinitely long- 
range and infinitely weak interactions" route to 
an explicit, analytic (mean- field- like) theory is well 
known jjlMlq]. What is new here is the connection be- 
tween nonintegrable interactions and the much-studied 
"nonextensive thermodynamics." This connection rests 
on an additional result, namely the demonstration that 
any ordering of the limits L, R — ► oo yields the same 
free energy. For the limit L — > oo first, the free energy 
is obtained directly from a Kac-potential treatment (at 
least for a < 0; the case a = is treated separately), 
while when the limit R — > oo is taken first, we obtain, 
as an intermediate step, a finite system with constant 
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interactions proportional to \/L d — i.e., without approxi- 
mation we obtain mean-field-like interactions! Finally, in 
the limit R oc L — > oo the independent prefactor multi- 
plying the pair interactions may be regarded instead as a 
power of L. This reproduces all scaling results of nonex- 
tensive thermodynamics of which we are aware, reducing 
the study of these systems to the application of stan- 
dard techniques with Boltzmann-Gibbs statistics. Fur- 
thermore, in the context of thermalized gravitational sys- 
tems |19]-£l| our results provide what we believe is a new, 
direct connection to a rigorous Kac-potential treatment. 
While a finite free-energy density is obtained in the 
R,L—too limit, the actual result depends explicitly on 
the cutoff junction. To be specific, consider a fluid of 
density p with the pair interaction 



oo r < a 

1 w{r/R) 



u(r) = < 



R d-T T T 

1 w(r/R) 
InR r^ 



r > a, <t < d 
r > a, t = d 



(1) 



where r = d + a is introduced to avoid confusion with 
the negative values of a, and where the cutoff function 
w(x) — taken to be isotropic for simplicity — decays at 
least as fast as l/x d ~ T+e for positive e, with w(0) fi- 
nite. The independent prefactors are chosen to preempt 
the divergence of the energy in the limit R — > oo. From 
this we obtain, in the limits L, R — * oo, the Helmholtz 
free-energy density 



f( P ,T) = C.E.{fo(p,T)-Ap 2 }, 



(2) 



where C.E.{. . .} represents the maximal convex envelope, 
/o (p, T) is the hard-core free energy (which is strictly pro- 
portional to T), and 



.4 



\S d $™ w(x)x d ~ 
\S d w{ty 



' l dx < r <d 

T = d 



(3) 



Here Sd = 2n d / 2 /T(d/2) is the surface area of a unit en- 
sphere. Up to a factor —ksT, A is just the second virial 
coefficient without the hard-core contribution. 

The free energy thus obtained indeed depends explic- 
itly on the cutoff function w(x), at least for r < d, but 
does not depend on the system shape. We obtain quan- 
titatively the same result for a lattice gas (hence also for 
Ising spin systems, cf. Ref. [|22f| . §11 E) with the substi- 
tution of the lattice hard-core free energy being the only 
modification. In this context, we recall that the first ap- 
plication of the Kac potential to spin systems is due to 
Baker (X6) . 

The dependence on details of the cutoff regulator has 
important implications when this solution is recast in the 
interpretation of "nonextensivity." For the non-periodic 
case of nonextensive thermodynamics, where the finite 
system size is used to regulate the energy, the bulk quanti- 
ties will depend on both boundary effects and the system 



shape, a point we have not found mentioned in previous 
studies. Furthermore, when periodic boundary condi- 
tions are employed and the interaction is cut off at some 
fraction of the system size, the bulk thermodynamics will 
depend on precisely which fraction is used, notwithstand- 
ing statements to the contrary || (cf. our discussion in 
Sec. IV below). The remainder of the paper is organized 
as follows. In Section [Il| we present our mathematical 
treatment of the various limits which lead to the results 
above. Next we present briefly the critical properties of 
these systems in Section III. Since we can interpret our 



results in the language of nonextensive thermodynamics, 
this connection is presented and discussed in Section IV, 



and various difficulties with nonextensive thermodynam- 
ics are brought to light, including the above-mentioned 
system-shape dependent "bulk" thermodynamics. Fi- 
nally, we conclude with a summary of our main results, 
and some remarks on the connection to work done in the 
area of gravitational systems. 



II. MATHEMATICAL TREATMENT 

A. Thermodynamic limit with finite range 

We begin by considering a fluid in d spatial dimen- 
sions with the pair- interaction potential given by (0), 
with finite R and in the thermodynamic limit. Then 
for < t < d our main results (||) and (||) for the limit 
R — > oo follow immediately from the rigorous treatment 
of Lebowitz and Penrose flS||. To see this, define 



<p(x) = w(x)/x T 



(4) 



in which case (Of) becomes u(r > a) = —R~ d (j)(r/R), the 
canonical Kac potential, with 4>{x) satisfying the nec- 
essary conditions for the proof given in p8[ , cf. Eqs. 
(1.21a)-(1.21c) in this reference. 

It is also possible and useful (for the r = d case) to 
understand this result from the Mayer cluster or virial 
expansion about a reference hard-core potential p2[ . 
Specifically, the Mayer function f(r) = e~P u ( r > — 1 (not 
to be confused with the free-energy density) may be de- 
composed as 



f(r) = [9(r - a) - 1] + 6{r - a)^- 011 ^ - 1] 



(5) 



with 0(x) the Heaviside step function. The first square- 
bracket term is the Mayer function for the hard-core po- 
tential and the second term accounts for the attractions. 
Each irreducible cluster of the virial expansion may be 
replaced with a sum of clusters in which the individual 
bonds are replaced, in turn, with each of the hard-core 
and attraction bonds. 

The hard-core bonds are independent of the cutoff R, 
so the sum of diagrams containing only these bonds yields 
the hard-core free energy (via Legendre transformation), 



which is unaffected by the limit R — > oo. In what fol- 
lows, we summarize the argument for why the attraction 
bond 2-cluster is the only other term which survives the 
R — > oo limit. 

The attraction bonds go to — 6{r — a)(3u{r) oc l/R d ~ T 
for large R (with r < d, for now). Each vertex that is free 
to integrate over space gives a factor J d d rw(r/R)/r T ~ 
R d ~ T . Consider first diagrams with only attraction 
bonds: the nth order ring diagram (each vertex having 
exactly two bonds) with n > 3 has n bonds and n — 1 free 
vertices, and so vanishes as l/R d ~ T for large R. All other 
n > 3 clusters have a higher ratio of bonds to vertices, 
and so these also vanish when all bonds are attractive. 

Replacing an attraction bond with a hard-core bond 
removes one factor of \/R d ~ T but also kills one of the 
free vertex integrals, i.e. the integration region is con- 
strained to the scale of the hard core. The net effect for 
ring diagrams is that they remain of the same order in 
R, and so vanish for n > 3. For more complicated irre- 
ducible clusters, a situation can arise where two vertices 
mutually constrained by hard-core bonds are also con- 
nected by an attraction bond. In this case, replacing the 
attraction bond with a hard-core bond does increase the 
order of the diagram by a factor R d ~ T . However, this 
situation can only occur for diagrams with higher powers 
of 1/R than the ring diagrams, and can never bring them 
up to order R°. Consequently, all diagrams with n > 3 
and at least one attraction bond vanish as R — > oo. 

The n = 2 case is distinct as it has one bond and one 
free integral, so it provides an R° contribution 



a 2 = hm - — 
H->oo 2 



,. -ps d 

lim 

R-^oo 2 



d*r Wm 



R d - 



w(x)x 



d—r—X 



dx 



a/R 



= -f3A 



(6) 



where A is given in (0). 

Thus the virial expansion reproduces (0), apart from 
the convex envelope. The failure of this otherwise exact 
method to reproduce the Maxwell construction is that 
the virial expansion relies on a homogeneous density and 
breaks down when this is not the case. Nevertheless, 
since it is exact whenever the density is homogeneous, 
the virial expansion supplemented with the second law 
(convexity) prescribes a unique free energy, and so can 
be regarded as providing the rigorous result. 

The utility of the virial-expansion method is that it 
applies to the borderline case r = d, where we cannot 
directly map to the results of Jig ]. The hard-core terms 
are unmodified, attraction bonds now go as 1/lni?, and 
free vertex integrals give 



d d rw(r/R)/r d -lni? 



(7) 



The sum of diagrams with n > 3 and at least one attrac- 
tion bond is 0(1/ In R) (provided there are no surviving 



resummations of l/(lni?) n terms). The n = 2 attraction 
bond diagram gives the only remaining contribution 



a 2 = hm -— 
R^,oo 2 



d d ru(r) 



lim — — — / w(x)x dx 



R^oc 2 In R 
= -±/3S d w{0) 



a/R 



(8) 



which again leads to (||) and (J3|). Interestingly, the bulk 
thermodynamics in this marginal case is not sensitive to 
the details of the cutoff function. 



B. Infinite range with finite system size 

Now consider the same fluid system, but taking R — ► 
oo with finite L. For clarity, we begin by considering a 
one-dimensional system with < t < 1, in which the 
periodic boundary conditions lead to the effective pair 
potential (r > a) 



u e g(r) = 



1 



w(\nL + r\/R) 



R l 



n=—oo 



\nL- 



(9) 



where r is understood to be less than L/2. As R becomes 
large, the direct (n = 0) interaction becomes negligible 
and an increasingly large number of terms contribute to 
the sum. By use of the expansions 



\n,L + r\- T = \nL\ 



and 



1 

nL 



-(t+1> 
2n 2 L 2 



(10) 



-r| |w£ ^" r * ) =w(\n\L/R) + sxgn(n)^w'(\n\L/R) 

+J^ w "(\n\L/R)±..., (11) 

where we assume w(x) to be analytic, we may rewrite 
Meff(^) for large R as 



u e g(r) 



R^\nL 

/ ' L_ tf_ 
+ \nR' R 2 



n L 
~R 



l + O 



(12) 



The convergence of this series is guaranteed by the shape 
of the cutoff function w and its derivatives. Note that 
the leading sum for large R is independent of the spatial 
separation r; the correction terms are down by a factor 
(L/R) 1 ~ T . In the limit of large R, we may, in turn, ex- 
press this sum as an integral 



lim u e ff(r-) = — - / (R/nL) T w(nL/R)dn (13) 

-R— »oo R 



LJo 



w(x)x T dx 



(14) 



For the borderline case r = 1, Eqs. (Q) and (O) have 
to be multiplied by (lni?) -1 . Furthermore Eq. (fi~3) now 
requires a lower integration limit n — 1 (or any constant) 
as a regulator. The pair interaction is then 



u e g(r) 



L\nR 



LlnR 



w(x)x 1 dx 



L/R 



[w(0)lnR + O(R )} 



(15) 



(16) 



In the limit R —* oo, we find the effective pair po- 
tential is again independent of the spatial separation, 
u eS (r) = -2L- 1 w(0). 

The generalization of this treatment to higher dimen- 
sionalities is straightforward. For a d-dimensional system 
of size L\ x L2 X ■ ■ • x L4, with periodic boundary condi- 
tions and < t < d, the effective pair potential (r > a) 
is 



UeffM 

where 



i? rf - 



• E • 

m — — 00 



E 



"■( 



f(ni,...,ra d ) 



f(ni,. 



,n d . 



7 (17) 



f(m, . . . ,rid) = |r + niLiXi + . . . + n d L d x d | (18) 



is the separation corresponding to the [n\, . . . , nj) peri- 
odic repeat, and the Xj are orthonormal vectors. In the 
limit of infinite interaction range i?, the direct interac- 
tion term (all m — 0) becomes negligible compared to the 
sum and we find the following generalization of Eq. (13): 



u e s(r) 



Sd 
V 



w(x)x dx , 



(19) 



where Sd has been introduced below Eq. (|3j) and V = 
Yli—i Li. We note that the strength of this constant in- 
teraction is inversely proportional to the volume, just as 
one expects for a mean- field- like system |23]. For t = d 
one has, in the limit R — > 00, u c s(r) = — SdV w(0). 
Since a constant pair interaction gives a configuration- 
independent energy density ^u gp 2 V, the energy and en- 
tropy contributions to the free energy decouple, and one 
obtains (0) directly in the thermodynamic limit. 

Instead of a "soft" cutoff provided by the (analytic) 
function w(x), one can also introduce a "hard" cut- 
off at a distance R, which is equivalent to taking, say, 
w(x) = 9(1 — x) in Eq. p). We start again with d = 1 
and t < d. The effective pair potential can be written as 



u e s(r) 



R 1 



R/L 

^ \n~L- 
=-R/L ' 



(20) 



where, for simplicity, we have taken R to be an integer 
multiple of L. For large R, the n = term becomes 
negli gibl e compared to the sum and by use of the expan- 
sion (|lC| ) we rewrite u c s(r) as 



2 (R 



-iR/L 

E 

71=1 



1 


l + O 


vW. 



(21) 



The 0(n~ 2 ) correction terms are again down by a fac- 
tor (R/L) T ~ l from to the leading sum [compare to 
(O)]. As before, we end up with a constant effec- 
tive interaction limfl^oo u e s(r) = —2L~ l /(l — r). For 
t = 1, the correction terms decay only like 1/lni? and 
lmiR^oo Meff(^) = —2/L. Since the case of general di- 
mensionality can be treated along the same lines, we only 
mention the resulting values for the constant pair poten- 
tial: lim^^oo u c g(r) = —SdV^ 1 /(d— r) for t < d and 
lmiR-joo M e ff(r) = — Sd/V for t — d. These results are all 
consistent with (||) with w(x) = 9(1 — x). 

Thus, we have shown that, in the limit of infinite in- 
teraction range, all pair interactions are identical even 
in a finite system, for any < r < d. For r > d on 
the other hand, the sum in Eqs. (g), ([17]) is convergent 
and hence the i?-dependent prefactor in u(r) should be 
omitted. The first term within the square brackets is 
then no longer negligible in the limit R — > 00, which 
leads to an effective interaction that depends both on r 
and on r and consequently to nontrivial critical behavior. 
Interestingly, the critical exponents retain their classical 
values until r > 3d/ 2, except for the correlation- length 
exponent v and the correlation- function exponent r\ [El. 
The critical temperature and other nonuniversal quanti- 
ties exhibit a nontrivial r dependence already for r > d, 
cf. Ref. fl24| and references therein. 



C. Lattice gas 

So far we have only discussed continuum fluid critical- 
ity. The generalization to a lattice gas (and, of course, to 
the Ising ferromagnet) results in the same free energy (0) 
with the appropriate lattice hard core fo(p). This may 
be seen most directly in the R — ► 00 first case, where 
the pair interactions go to a constant. This derivation 
applies both for particles in a continuum or on a lattice 
(indeed, no specification was made), and so the constant 
interactions are obtained in both cases. The final step of 
constructing the free energy from the decoupled energy 
and entropy reveals that the lattice hard-core free energy 
is the appropriate one to use in (0). 

For the L — > 00 first limit, one must use a lattice gen- 
eralization of the Mayer expansion. Such an expansion 
would give for non-interacting particles the lattice hard- 
core free energy. The lattice sums for clusters with at- 
traction bonds can be taken to integrals in the large R 
limit, and so the rest of the continuum derivation applies. 



D. Thermodynamic limit with a system-size 
dependent range 

The final case we consider is that of R ex L — > oo, the 
case most directly applicable to "nonextensive thermody- 
namics." Consider a finite system with periodic bound- 
ary conditions (with all dimensions Li oc L) and direct 
pair interactions given by (fil) for some finite R. The pe- 
riodic boundary conditions give rise to the effective pair 
potential (JIT]) . The virial expansion for the pressure may 
then be obtained from a finite-volume cluster expansion 
with the Mayer function / = cxp(— (3u e s) — 1 ~ —0u e s 
for large R. Since the effective interaction is a sum of pair 
interactions, a simple u c g cluster such as the irreducible 
3-cluster shown in Fig. |l|(a) decomposes into a large num- 
ber of pair-interaction bond clusters [cf. Fig. 0(b)]. Note 
that, in order to avoid overcounting, one end of every 
bond must remain in the original finite volume, but the 
other end may be taken to lie in any of the periodic re- 
peats (and these still contribute to the virial expansion 
as they are parts of an irreducible cluster) . 




(a) (b) 

FIG. 1. (a) The 3-cluster formed with the effective inter- 
action u c s and (b) one contribution to the decomposition of 
this cluster into pair interactions u(r). 

In spite of these complications, we remark that all 
hard-core bonds are the same as in the previous (L — > oo 
first) case, since they only appear in the direct interac- 
tion of the original pair. Hence the hard-core bonds sum 
to the hard-core free energy regardless of whether L and 
R are simultaneously taken to infinity, or L first. 

Furthermore, all n-clusters with n > 3 and at least 
one attraction bond can be shown to still vanish as 
R,L—> co. First, all bonds still carry a factor of l/R d ~ T 
for large R, as before. Free vertices integrate to 



d d 



r r T w(r/R) 



R 



d-7 



(L/R) d 



d d xx~ T w(x) (22) 



for the direct interaction of the original pair, which goes 
as R d ~ T times a finite factor for large R. For a pair in- 
teraction involving a neighboring replica, the vertex in- 
tegration is the same as above with a shift in the argu- 
ment of w(f/R), where f [see Eq. (18)] goes roughly as 
r + cL. The resulting vertex integration will also scale as 
R d ~ T , provided r + cL < R. Since f increases for each 
increasingly remote replica, the cutoff function w(x) will 
ensure that only a number 0((R/L) d ) of such bonds will 



contribute. Hence free vertex integration with effective 
interaction bonds, while considerably more complicated, 
still results in a fact or R d ~ T at large R. Previous argu- 
ments from Section II A then apply, and so these terms 
all vanish when R,L—* oo. 

The remaining 2-cluster integral may be written as 



«2 



f 
lim - / d ru e g(r) 
R-*oo Z /_^„ 



w(f/R) 



(23) 
(24) 

. (25) 



where the sums cover all replicas (and we have omit- 
ted writing the hard-core condition for clarity). In going 
from (E4|) to (J25J) the integrand variable changed from f 
to r, consistent with the definition of f. The remaining 
integrals piece together a single volume integral over all 
space, so 




« 2 



lim 



-/3 



2R d - 



d d rw{r/R)/r 



-0A, 



(26) 



following (|6j). Combining this with the hard-core contri- 
bution gives the same free energy (0) in the R, L — ► oo 
limit as was found previously for L or R going to infinity 
first. The significance of this simultaneous limit towards 
nonextensive thermodynamics will be discussed in Sec- 
tion Ev. 



III. CRITICAL BEHAVIOR 

The critical properties of nonintegrable systems are 
readily found from the Hclmholtz free-energy density (0) 
via standard procedures. For example, the critical den- 
sity can be obtained from d 3 f(p,T)/dp 3 \ Pc — 0, which 
reduces to the temperature (and attraction) independent 
condition d 3 fo(p,T)/dp 3 \ Pc = 0. The critical temper- 
ature is then found from d 2 f(p,T c )/dp 2 \ Pc — which 
gives 



k R T r = 2A 



d 2 (Pfo) 
dp 2 



(27) 



For a hypercubic lattice with lattice constant a, the 
lattice-gas hard-core free energy is P/q = p In p+ (a~ d — 
p) ln(l — a d p). This results in the critical values 



lattice gas: a d p c = 1/2 k B T c = A/(2a d 



(28) 



Although the free energy for continuum hard spheres (di- 
ameter a) is not known exactly, a very good approxi- 
mation in three dimensions is nevertheless given by the 
Carnahan-Starling (CS) expression 



/3/o° b (p) = /3/ M (p) 



p 2 v {4 - 3pv Q ) 
(1 - pvo) 2 



(29) 



with wo = 7ra 3 /6 the volume of the hard sphere and / Id 
the ideal gas free energy density. Using this for f yields 
a quintic equation for p c with a unique positive root ]2q] 
and the critical values 



CS: 



a 3 p c 



0.249129 k B T c ~ 0.180155A/a 3 . (30) 



The susceptibility %, defined as the ratio of the isother- 
mal compressibility to the ideal gas compressibility, 



X -\p,T) = p 



d 2 ((3f) 
dp 2 



(31) 



exhibits near criticality the classical divergence \ ~ C+/£ 
for p — p c and positive reduced temperature t = (T — 
T c )/T c . From @ we find 



C*4 



PC 



d 2 (Pfo 



dp 2 



(32) 



which results in C+ G = 1/2 for the lattice gas and 
C9 S ~ 0.361569 for the d = 3 hard-sphere fluid, ap- 
proximated by the CS expression. 

Finally, the order parameter for T < T c is given by 



a d \p-p c \ = BV~t 
for t sufficiently small, with 



B 



„, , 18[d 2 f /dp 2 ] Pc 

[d 4 h/d P % c 



1/2 



(33) 



(34) 



This evaluates to B LG = 3/2 and B 



? OS 



1.13459. 



IV. NONEXTENSIVE THERMODYNAMICS 

As indicated in the Introduction, the main motivation 
for this work stems from the considerable attention sys- 
tems with nonintegrable interactions have received in the 
context of "nonextensive thermodynamics." An essential 
aspect of these studies of nonextensivity is the use of the 
system size as the regulator for the energy. Furthermore, 
the interactions are not scaled by a negative power of the 
system size but left with strength of order unity. 

Our cutoff interaction with range R can be interpreted 
directly in terms of this nonextensive thermodynamics for 
R/L constant, which reproduces the system-size regula- 
tor. The negative power of L multiplying our pair inter- 
action is eliminated by scaling the temperature according 



to T -> L d - T T in the Boltzmann factor (T -> Tln£ for 
t — d). The scaling of this system-size dependent "tem- 
perature" (and consequently the free energy) matches ex- 
actly the conjectures of nonextensive thermodynamics, 
thus we have derived rigorously the primary conclusions 
of Refs. |g||[ll],[12J using only standard methods. We 
note that, despite what has been suggested in Ref. pq , 
neither the explicit free energy nor exact results for the 
nature of the criticality were obtained in Refs. fuyj- 

Next we want to discuss some examples from the recent 
literature and point out a few problems attached to the 
interpretation of nonextensive thermodynamics. First, a 
pervasive notational problem in the nonextensive ther- 
modynamics literature is the use of "long-range interac- 
tions" to mean "nonintegrable interactions." The former 
term already has a standard meaning within the con- 
siderably more important class of integrable interactions 
p7| , p8| Jl|. The converse problem also exists, where all 
integrable interactions (including the true long-range in- 
teractions) are termed "short-range" (see, e.g., Ref. p9[). 

Many papers addressing nonintegrable interactions 
have relied on numerical simulations of finite systems. In 
general, the required regularization of the interactions is 
carried out either by imposing free boundary conditions 
(which leads to an inhomogeneous system) or by adopt- 
ing periodic boundary conditions and cutting off the in- 
teraction at half the system size, being the maximum 
separation between the particles. For example, Curilcf 
and Tsallis || have performed molecular-dynamics sim- 
ulations of fluids in d = 2,3 dimensions with Lennard- 
Jones-like interactions, with an attractive tail decaying 
like r~ T and 1 < r < 2d. It is explicitly stated (Ref. M, 
p. 271) that, in the thermodynamic limit, no physical 
consequences emerge from the ( "computationally conve- 
nient" ) adoption of a cutoff at half the size of the simu- 
lation box. Our exact solution presented in Sec. |ll| shows 
that this is not correct for r < d, but that rather dif- 
ferent bulk thermodynamics emerges for different cutoff 
distances and cutoff functions. 

In Ref. |g9[, rings of magnetic particles in a col- 
loidal suspension have been studied numerically, where 
the nonmagnetic part of the interactions has the above- 
mentioned generalized Lennard-Jones form. The main 
results are heuristically interpreted in terms of "nonex- 
tensive thermodynamics" by observing that the size de- 
pendence of the total energy of the rings can be well 
described by a scaling law obtained from the integrated 
interaction (which is essentially a mean-ficld-like approx- 
imation). This scaling law is just what is also found 
within the "g-generalized thermodynamics" of Ref. |j|, 
commonly referred to as Tsallis g-statistics, and conse- 
quently a r dependence (r = d + a) is proposed for the 
so-called "nonextensivity parameter" q appearing in this 
formalism: q = \ (corresponding to Boltzmann-Gibbs 
statistics) for r > d and q — 2 — d/r for t < d. Also 
in Ref. 1 30 1 it has been conjectured that q is a d- and 
r-dependent quantity for r < d. It appears that these 



conclusions have been since revised 0,^1) , and these sys- 
tems are now classified as "weakly violating" Boltzmann- 
Gibbs statistics, meaning that q = 1, whereas several 
thermodynamic quantities lose their extensivity. Never- 
theless, recent studies have continued to explore q =/= 1 
values for r < d, citing an alleged natural connection 
between Tsallis statistics and nonintegrable interactions 
|p2| . Our results show that the same energy scaling for 
t < d may be obtained with the conventional value q = I. 
Inspired by the system-size dependence of the energy 
found in Ref. J29|j , an Ising model with interactions de- 
caying as r~ T (r > 0) and free boundary conditions has 
been analyzed by mean- field methods in Ref. || . On the 
basis of the resulting values of the critical temperature for 
t = and r = d, it was then conjectured that the mean- 
field prediction for the critical temperature might hold 
for all < t < d. Our exact result now demonstrates 
that this is indeed true, but reveals in addition an awk- 
ward consequence of the boundary conditions adopted 
in Ref. M: since the thermodynamic limit explicitly de- 
pends on the choice of the cutoff function, an inhomoge- 
neous system with an inhomogeneous cutoff will lead to 
bulk thermodynamics that depends on the shape of the 
system, which to our judgment constitutes an undesirable 
feature of the nonextensivity formulation. We note that 
already in an early study of the Ising model with long- 
range interactions Q a system-shape dependence of the 
thermodynamic properties has been observed for the case 
of conditionally convergent lattice sums with r = d (such 
as dipolar forces). Also the exactness of mean- field the- 
ory for shape-independent forces in the limit r — ► d+ has 
been obtained in the same reference, essentially from an 
observation similar to ours for general < r < d, namely 
that one divergent term dominates all other terms in the 
lattice sum. However, for the case r < d, we have been 
unable to find in the literature any mention of the system- 
shape dependent thermodynamics that must result for all 
inhomogeneous systems with nonintegrable interactions 
and a system-size dependent cutoff. 



V. SUMMARY AND CONCLUSIONS 



In summary, we have shown that (i-dimensional, peri- 
odic systems with nonintegrable, algebraically decaying 
interactions, i.e., interactions of the form u(r) ~ — l/r r , 
with < t < d, are exactly described by mean-field the- 
ory, upon introduction of a cutoff R in the interaction 
range and the proper independent rescaling of the in- 
teraction strength. This proof holds for either order of 
limits R — > oo and L — > oo (where L is the linear system 
size), including the simultaneous limit, and the resulting 
free energy depends on the details of the cutoff. 

Our study employs Boltzmann-Gibbs statistics and 
pertains directly to nonextensive thermodynamics, pro- 
viding explicit, exact results for the thermodynamics and 
critical behavior. In doing so, we show that nonintegrable 



interactions do not require |3C],|32| the application of gen- 
eralized g-statistics. Furthermore, the explicit regulator 
dependence — cutoff length, cutoff shape, and even sys- 
tem shape for inhomogeneous systems — is demonstrated, 
a topic which has been mostly neglected in nonextensiv- 
ity studies. 

On an intuitive level, our findi ngs for the case of finite 
systems with infinite R (Section [II B ) result from the di- 
vergence of the lattice sums over the periodic copies of the 
system under consideration; these divergent sums then 
dominate the direct pair interaction, and are, to leading 
order, independent of the spatial separation between par- 
ticles. A suitable normalization is indispensable for the 
existence of the thermodynamic limit and — as we have 
pointed out — a regulator depending on the number of in- 
teractions rather than on the system size emerges as the 
natural choice. The resulting effective pair interaction is 
then independent of the spatial separation, for < r < d. 

We have concentrated on both fluids and lattice gases. 
However, since our method for the R — > oo first case 
demonstrated how the effective interaction becomes in- 
dependent of spatial separation, our results immediately 
carry over to large classes of other systems as well (obvi- 
ously, the critical properties obtained in Sec. [II explicitly 



refer to systems with a one-component order parameter). 
These include general O(n) models (XY, Heisenberg, . . .) 
and Potts models. In this context we note that in ear- 
lier work the exactness of mean-field theory in the limit 
t — > d+ has been found not only for the Ising model 
(n = 1) |33|] but also for the Husimi- Temper ley mean 
spherical model (n — ► oo) 34 1. 

Our explicit result of an analytic free energy, gener- 
alized to the systems mentioned above, explains a num- 
ber of numerical results obtained for systems with non- 
integrable interactions. This includes the molecular- 
dynamics simulations of systems with a generalized 
Lennard- Jones potential discussed in Sec. IV (3o> Morrtc 
Carlo simulations of one-dimensional Ising [|9| , Potts [jllj , 
and XY |IJ] models with r < 1, and the scaling proper- 
ties found in Refs. @|l(p|,ll|] ■ In Ref. pi, a subleading 



term in the spin-spin correlation function of the nonin- 
tegrable Ising chain was considered and on the basis of 
Monte Carlo simulations of finite systems with t = 0.50 
and t — 0.75 it was concluded that these correlations are 
correctly described by mean-field theory. 

Finally, we remark on the connection of these results 
to gravitational systems, where d = 3 and r = 1. Since 
the masses of the particles and the gravitational cou- 
pling are presumably fixed, we are not at liberty to scale 
to infinitely weak interactions. However, following the 
"nonextensive thermodynamics" formulations, we can 
consider R oc L for a finite system and regard the pref- 
actor in the pair interactions as belonging to a rescaled 
temperature T — > L d ~ T T = L 2 T. At the same time, we 
note that such a divergent temperature appears of lim- 
ited practical use. 

The existence of a phase transition in these sys- 
tems has been studied for fermionic particles in the 



context of Thomas-Fermi theory |1^,^(J. Most stud- 
ies have considered gravitationally interacting particles 
in a non-periodic, finite-sized system, concentrating on 
the asymptotic, large-iV limit. Here they obtain non- 
extensivity in the particle number N, with the energy 
per particle growing as iV 4 / 3 (in contrast to TV for in- 
tegrable interactions). In Ref. [£1[ the energy densities 
were suitably rescaled to enable the infinite- volume limit, 
which presumably comes closest to what has been pre- 
sented in the current work. 

However, these studies differ from ours due to the 
fcrmionic character of the particles, which is used to reg- 
ulate the short-distance behavior. Classical gravitational 
systems have previously been studied as well p7J , |38| | , 
again in a non-periodic, finite-sized box, with various 
forms of short-distance regulators. In this case the en- 
ergy per particle scales as N, essentially because the pair 
interaction, while decaying as 1/r, has a minimum value 
for all particles that is proportional to 1/L, due to the 
fixed system size. In contrast, our thermodynamic limit, 
L — > oo with fixed particle density, combined with the 
long-distance power-law tails of the pair potential, im- 
plies an energy per particle growing as N 2 / 3 (or N 1 ~ T / d 
for general r < d). Finally, we mention Ref. g^| for a 
recent review of some interesting features of thermalized 
gravitational systems. 
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